Hidden order in one dimensional Bose insulators 
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We investigate the phase diagram of spinless bosons with long range (oc 1/r 3 ) repulsive inter- 
actions, relevant to ultracold polarized atoms or molecules, using DMRG. Between the two con- 
ventional insulating phases, the Mott and density wave phases, we find a new phase possessing 
hidden order revealed by non local string correlations analogous to those characterizing the Haldane 
gapped phase of integer spin chains. We develop a mean field theory that describes the low energy 
excitations in all three insulating phases. This is used to calculate the absorption spectrum due to 
oscillatory lattice modulation. We predict a sharp resonance in the spectrum due to a collective 
excitation of the new phase that would provide clear evidence for the existence of this phase. 



Introduction - Systems of ultracold atoms in optical 
lattices hold significant promise in the study of corre- 
lated quantum matter. Experimental realization of the 
Bose-Hubbard model has been an important step in this 
direction, not least because it facilitated the first obser- 
vation of a quantum insulating state of bosons^. Lo- 
calization of bosons in the Mott insulator is driven by 
strong on site interactions leading to an extremely sim- 
ple state, well described by a site-factorizablc mean field 
wave- function. An interesting question is what other 
phases, perhaps with non trivial structure, may be stabi- 
lized by longer range interactions 2 ^^^. This issue has 
been given added urgency by recent advances in trapping 
and cooling of atomsi and molecules^ with large dipole 
moments. 

Ultracold dipolar bosons, in one dimensional optical 
lattices can be described by the hamiltonian 
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where we assumed that the dipoles arc polarized by an 
external field perpendicular to the lattice. The three pa- 
rameters appearing in the hamiltonian can be tuned in- 
dependently in experiments 2 -. The conventional phases of 
this system at integer filling are known from mean field 
studies 2 -^. They include the Mott insulator (MI) at large 
U, a density wave (DW) for large V, and a supcrfluid 
(SF) for large t. 

In this letter we use the Density Matrix Rcnormaliza- 
tion Group (DMRG) method 9 - to show that a new gapped 
insulating phase with non trivial structure obtains in a 
wide parameter regime between the two conventional in- 
sulators. The new phase is separated from the conven- 
tional insulators by lines of second order transitions and 
is found to possess particularly subtle order that is re- 
vealed only by highly non local string correlation func- 
tions. Correlations of similar nature exist in the Haldane 
gapped^ phase of quantum spin-1 chain a 11 ' 12 . We shall 
therefore term the new phase as the Haldane Bose insu- 
lator (HI). The analogy can be made more explicit by 
truncating the Hilbert space of the Bose system to three 



occupation states per site (for example, for a system with 
ft = 1 particles per site, we keep only the occupation 
states n = 0,1,2 for every site). This defines an effec- 
tive spin-1 model with S* = rii — n and n the average 
filling. In this space the DW phase corresponds to anti- 
ferromagnetic ordering of the pscudospins in the z direc- 
tion. The MI ground state, on the other hand, includes 
a large amplitude of the state with Sf = on every site 
and small admixture of states containing tightly bound 
particle- hole fluctuations (S z = ±1 on nearby sites). The 
Haldane phase may also contain mostly sites with S z = 
at large U, but the ordering of the fluctuations is unusual. 
The string correlations imply that particle and hole fluc- 
tuations appear in an alternating order along the chain 
separated by strings of zeros of arbitrary lengt h 11 ' 12 . The 
truncation to three states is justified only at large U when 
fluctuations in site occupancy are strongly suppressed. 
However, from the numerical results we conclude that 
the long ranged string order of the actual bosons sur- 
vives fluctuations in the occupancy beyond the effective 
spin-1 description. 

After a description of the DMRG results, we will show 
how the new phase can be detected experimentally. The 
challenge lies in the fact that standard experimental 
probes couple to local observables and are blind to the 
highly non local string correlations. We shall therefore 
look for signatures of the new phase in the excitation 
spectrum rather than the ground state properties by con- 
sidering the response to lattice modulation^. We de- 
velop an approximation scheme that enables us to calcu- 
late the response functions in all three insulating phases. 
We predict a sharp resonance in the absorption due to a 
neutral collective mode (S z — in the pseudospin termi- 
nology), that is special to the Haldane insulator and can 
serve to identify this phase. 

Numerical Results - We investigate ground state and 
lowest excitations of the hamiltonian |T|) in a space of 
the parameters U and V using the DMRG algorithm 9 . 
Throughout we consider a filling of n = 1 particles per 
site in the ground state. The maximal length of the sys- 
tem is N = 256 sites, and up to M = 200 states are kept 
per block. In most of the parameter regime we found that 
the results do not change significantly when the cutoff in 
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boson occupation number per site is increased beyond 
4. The l/?' 3 interactions was included up to the next 
nearest neighbor range. Open boundary conditions were 
used, and opposite external chemical potentials were ap- 
plied on the first and last site in order to lift the ground 
state degeneracy in the HI and DW phases. 

The phase diagram of ([I]) in the (U, V) plane is shown 
in Fig. [1] The nature of the phases in the DMRG simula- 
tion was elucidated by a direct calculation of the ground 
state correlation functions: 

Rs F (\i-j\) = (btbi) (2) 

^Dw(H-j|) - (-I) 1 * -41 (SmSnj) (3) 

String (|i-j|) = (Sm^i^Hn^ (4) 

Here Sn.i = rii — n. The supcrfluid phase is character- 
ized by a power law decay of Rsf oc \i — j| -1 ^ 2 ^ with 
Luttinger parameter K > 2. In the MI phase, all the 
correlation functions decay exponentially to zero. In the 
DW phase Rdw (|* ~ j\) approaches a constant at long 
distances (note that also i? s tring(|i — j\) — ► const. ^ 0), 
and the lattice translation symmetry is spontaneously 
broken. In analogy to the S = 1 XXZ chainilHii^, 
we expect the appearance of another phase (the Haldane 
phase) between the MI and DW phases. This phase is 
characterized by i? s tring (|* — j\) const. ^ 0, while the 
DW correlations decay exponentially Unlike the den- 
sity wave phase, this phase does not break the lattice 
translation symmetry It does, however, break a hidden 
Zi symmetry related to the string order paramete r 12 ' 16 . 
Fig. [21a) presents an example of how the phase bound- 
aries were determined: we show the string and DW order 
parameters (defined as the square root of the asymptotic 
values of the corresponding correlation functions) as a 
function of V along the line U — 6t. The phase tran- 
sitions from MI to HI and from HI to DW are clearly 
visible and seem to be of second order. 

The phase diagram does not change qualitatively if we 
keep only the nearest neighbor interactions in JT]). How- 
ever, further range interactions act to frustrate the DW 
order and thereby widen the domain of the HI phase. We 
note that previous DMRG studies of the Bosc-Hubbard 
model with nearest neighbor interaction [17] did not look 
for the string correlations and therefore did not find the 
subtle HI phase. 

In addition to the ground state, the energies of the 
first few excited states were calculated. The gap to the 
first excited state with the same number of particles as 
the ground state, A = Eg^ =Q — £$° =0 , was calculated 
by targeting also the first excited state in the DMRG 
calculation. Here 5n is the number of particles relative to 
a state with exactly n = 1 particles per site. The charge 
gap of the system Ai = E}°J =1 + E% ) J = _ 1 - 2E$ =0 was 
calculated by targeting the ground states of the Sn = ±1 
sectors. 

The gaps Ao, Ai along the line U = 6t in the (U, V) 
plane are shown in Fig. [U^b). At the transition point 
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FIG. 1: Phase diagram of the hamiltonian ([TJ in the {U,V) 
plane, obtained from DMRG. The phases that appear in the 
diagram are: Superfluid (SF,o), Mott insulator (MI,D), den- 
sity wave (DW,x) and Haldane insulator (HI, A). INSET: 
density wave (dashed line) and string (solid line) correlations 
at a particular (U,V) point (see (|3l4fl h This point is in the 
HI phase, as can be seen from the fact that string correlations 
do not decay, whereas DW correlations decay rapidly. 



between the MI and HI phases, both Ao and Ai vanish, 
while at the transition between the HI and DW phases 
only Ao vanishes. This indicates that both transitions are 
second order, but the nature of the critical excitations at 
the transition is different. 

In the MI phase, the lowest excitations are particle 
and hole excitations. This is clear from the fact that 
Ai = Ao, i.e., the first excitation with 5n = is an 
unbound particle-hole pair. The gaps to these particle 
and hole excitations vanish at the transition to the HI. 
The Haldane phase displays another low-energy excita- 
tion. In a certain region of the phase diagram Ao goes 
below the charge gap indicating the presence of a gen- 
uine neutral mode with 8n = 0. In Fig. HJb) the crossing 
point between the two states is clearly seen as a cusp 
in the Ao curve near the middle of the HI phase region. 
At the transition to the DW phase, only the neutral gap 
vanishes. The presence of the neutral mode seems to be 
characteristic of the Haldane phase, and can be used to 
detect it, as will be discussed below. 

Experimental detection - In the numerical simulation 
we were able to establish the presence of the string or- 
dered insulating phase HI by measuring the non local 
ground state correlations (J4| directly. By contrast, exper- 
imental probes naturally couple to local operators, such 
as the charge density. Our detection strategy will focus 
on probes of the excitation spectrum, which may exhibit 
distinct (albeit indirect) signatures of the HI phase. 

Let us consider in some detail the response to paramet- 
ric excitation of the optical lattice. This technique was 
used successfully to probe the excitations in both the MI 
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FIG. 2: (a): DW order parameter (A) and string order pa- 
rameter (□) as a function of V along the line U — 6t. (b): 
Gaps to neutral and charged excitations, Ao (o) and Ai (o) 
respectively (see text) along the same line. 



and SF phases^. The idea is to apply a periodic modula- 
tion of the lattice intensity. Within the lowest Bloch band 
this perturbation simply corresponds to uniform modula- 
tion of the hopping matrix element in ([1]) . The perturbed 
hamiltonian can be written as H + h cos(ujt)T, where T is 
the kinetic energy operator T = ^ b\bi+i+H.c.. Within 
linear response, the absorption rate is 
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We shall calculate I(u>) within the effective spin-1 hamil- 
tonian 



H = ^E(^r + i+H.c.) 
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using the corresponding perturbation operator T = 
J"V S^ +1 +H.c. The mapping to a spin model amounts 
to projecting ((T|) on the subspacc including only three oc- 
cupation states per site (5rii = 0, ±1). This is justified 
in the insulating phases at large U, where multiple par- 
ticle or hole occupancy is suppressed. The U term in ([6|) 
corresponds to on site interaction, V to nearest neighbor 
interactions, and J to boson hopping. The projection 
of the hoppnig term in (JTJ) gives another contribution to 
the spin model, that breaks the particle-hole symmetry. 
Omitting this term, as we have done in ([6]) , does not lead 
to qualitative changes in the phase diagram. The advan- 
tage in using the model (JB]) is that it is amenable to a 
mean field theory, developed by Kennedy and Tasakii^, 
that captures all three insulating phases. 

Kennedy and Tasaki introduced a non local unitary op- 
erator that transforms the string correlation ^ to con- 
ventional spin correlations {S^Si-}, which admit a local 
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FIG. 3: Absorption spectrum ([5]l calculated from the mean 
field theory along the dotted line in the inset (mean field phase 
diagram of (JTJl ) . The dotted line in the spectra marks the 
position of a quasi-particle peak seen only in the HI phase 
(b). The dashed and dash-dotted lines mark the edges of two 
particle continua. 



mean field treatment. At the same time the hamiltonian 
flS]) assumes a rather unusual, but nevertheless local form: 

h = ./^.v;.s;. , s'!cx V a-s: • /-.v;. ; ;.s7. ; 



(7) 



The U(l) symmetry of the original hamiltonian seems 
to have been broken down to x Z2. In fact the full 
U(l) symmetry is preserved but is generated by a highly 
non-local operator in the new representation. 

The mean field theory of (J7J) consists of finding the 
best site-factorizable wave-function for this hamilto- 
nian. The solutions are of the general formic | } = 
]li(cos0 I 0) i ± sin# I ± 1 } i ). The resulting phase dia- 
gram is shown in the inset of Fig. [3J In the Haldane phase 
< 9 < tt/2 and the ground state is four fold degenerate, 
which reflects broken Z2 x Zi symmetry. 9 is given explic- 
itly by the expression cos 2 6 = (U - 2V + 4J)/(8 J - 2V). 
The DW phase is characterized by sin# = 1, i.e. dou- 
bly degenerate ground state. Finally in the MI there is 
a non degenerate ground state with cos6> = 1. The fact 
that the supcrfluid phase is pushed to negative V is an 
artifact of the projection to three occupation states. 

The collective excitations are found, following a stan- 
dard scheme^, by quantizing the small fluctuations 
around the mean field minima. We obtain an effective 
hamiltonian of the form H e ff = w a kfy Jfiak for the 
two collective modes (3 a k (a = 1,2). Likewise we ex- 
pand the perturbation T in terms of the collective mode 
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operators. Finally we compute I (to) using ([5]). 

Fig. [3J displays the results at different points along 
a line in the (U, V) plane that cuts through the three 
phases. As expected we see the gap in the spectrum clos- 
ing at the two quantum phase transitions. The most dis- 
tinct signature of the Haldane phase is a delta-function 
peak in the absorption spectrum corresponding to cre- 
ation of a single quasi-particle. This resonance can serve 
as a "smoking gun" experimental evidence of the HI 
phase. The fact that this quasi-particle couples to T 
which does not change particle number suggests that it 
must be a neutral excitation (quantum number Sn = 0). 
Similarly we infer that the other quasi-particle, which 
does not couple to T, is a charged mode. The other two 
peaks seen in the spectra are the bottom edges of the two 
particle continua corresponding to a pair of "charged" 
and a pair of neutral excitations respectively. 

While the mean field treatment is approximate and 
cannot be expected to capture details such as the loca- 
tion of phase boundaries, it agrees with the DMRG re- 
sults on the qualitative features. In both calculations the 
signature neutral mode is below the particle hole contin- 
uum in a wide region inside the HI phase (see Fig. [2](b) 
and Fig. [3J gray region in the inset). The mean field 
theory also captures the fact that only the neutral mode 
becomes critical at the transition from the HI to the DW 
state. Because of the lattice symmetry breaking in the 
DW phase, we expect a neutral excitation with momen- 
tum k = 7r to drop to zero energy at the transition to 
this phase. This mode is adiabatically connected to the 
S* = magnon of the spin-1 Heisenbcrg chain. However 
it cannot be identified with the delta function peak in Fig. 
|3j because the operator T couples only to k — excita- 
tions. The delta function peak is therefore interpreted as 
a k = bound state of two such k = ir excitations. A 
single fc = 7r excitation will show up as a sharp peak in 
Bragg spectroscopy^, in which the lattice is modulated 
with another laser of wave-vector k = n. 

We now discuss the possibility of realizing the HI phase 



experimentally in a system of ultra-cold atoms in a one 
dimensional optical lattice. From Fig. Q] we see that 
this requires V ~ U. This condition can be expressed 
in terms of the elastic scattering length, a s , and the 
effective scattering length associated with a dipole d, 
cid = —md 2 /4:irh 2 . For a typical optical lattice poten- 
tial, with wavelength k^ and effective depth of the or- 
der of 10E R (E R = (hk L ) 2 /2m), we estimate^ U w 
l2(kLa s )Eji. The nearest neighbor interaction energy is 
V = d 2 /r 3 w d 2 kl/ir 3 and therefore the condition V ~ U 
is equivalent to \aa\ ~ 30a s . Since the typical scattering 
length of atoms is of the order of tens of angstroms, we 
need an effective dipole scattering length of order 1000 A. 
Dipolcs of this magnitude can be obtained with ultracold 
dipolar molecules^. Alternatively it might be possible to 
use Cr atoms with w QA if the other energy scales t 
and U are reduced proportionally. Of course the temper- 
ature would have to be lowered to the same scale, below 
the gap A ~ t ~ V seen in Fig. Ufa). This is achieved 
automatically by adiabatic cooling as the optical lattice 
potential is increased slowly to the desired value. The 
fact that the HI phase is incompressible implies that in 
the presence of a harmonic confinement it will form a 
plateau structure, similar to the usual Mott insulator. 

Discussion and conclusions - We have shown that ad- 
dition of further range interactions to the Bosc-Hubbard 
model can give rise to a new insulating phase with non lo- 
cal string order, analogous to the Haldane gapped phase 
of integer spin chains. The new phase can be realized in 
systems of ultracold dipolar atoms or molecules in opti- 
cal lattices. We predicted unique signatures of the new 
phase in the parametric excitation spectra. These include 
low energy critical excitations that appear near the phase 
transitions to the conventional phases, and a sharp reso- 
nance in the response seen only in the new phase. 
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